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Abstract 


Design  and  packaging  of  vehicle  interiors  and  cockpits  had  become  a  science  on  itself, 
particularly  in  recent  years  where  safety  is  paramount.  Significant  experimental  efforts 
have  been  made  to  study  driver  reach  and  barriers  as  they  have  direct  effect  on 
performance  and  safety.  This  paper  presents  a  rigorous  formulation  for  addressing  the 
reach  envelope  and  barriers  therein  of  a  3 -point  restrained  driver  and  compares  with  a 
lap-restrained  driver.  The  formulation  is  based  on  a  kinematic  model  of  the  driver, 
which  has  the  upper  body  and  arm  characterized  as  a  7  degree  of  freedom  (DOF)  for 
unrestrained  and  4DOF  for  3-point  restrained.  Those  kinematic  equations  are  further 
developed  to  address  crossability  analysis.  Visualization  of  such  barriers  and  their 
crossability  results  within  the  reach  envelope  provide  significant  insight  into  driver 
performance  and  reach  zones. 


Keywords:  SAE,  driver  reach,  ergonomics,  workspace,  vehicle  interior  design 
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Introduction 


Industrial  tasks  often  involve  repeated  operator  reach  to  various  tools,  parts  and  controls 
placed  around  a  seated  workstation.  Examples  include  aircraft  cockpit  design,  vehicle 
interior  design  and  manufacturing  workstations.  All  controls  and  other  elements  in  the 
workstation  requiring  frequent  manual  operation  should  be  placed  in  an  area  that  focuses 
on  optimizing  the  human  interface  for  the  target  population,  minimizing  occlusion  and 
reach,  while  providing  maximum  leverage  for  operation.  The  driver  reach  capability  is  a 
basis  for  developing  a  design  aid  in  the  form  of  reach  contours  that  can  be  used  easily  by 
designers  during  the  conceptual  phase  of  a  new  vehicle  program.  While  a  database  of 
reach  barriers  can  be  tabulated  into  a  large  database,  we  envision  the  implementation  of 
the  methodology  into  a  digital  human  and  simulation  program  such  that  reach  barriers  are 
calculated  and  visualized  when  needed  during  the  design  process. 

The  tool  used  for  representing  driver  reach  capabilities  is  the  model  developed  in  SAE 
Recommended  Practice  J287  by  Hammond  and  Roe  (1972).  The  SAE  surfaces  are 
parameterized  by  a  packaging  factor  that  combines  vehicle  interior  dimensions  such  as 
seat  height,  fore-aft  and  vertical  steering  wheel  position,  and  steering  wheel  diameter, 
into  a  single  “G”  score.  This  model  is  based  on  measurements  taken  of  the  actual  reach 
capabilities  of  a  representative  sample  of  drivers  and  accounts  for  the  design  variation  in 
workspace  geometry  of  a  wide  range  of  vehicles  in  term  of  seat  location,  steering  wheel 
position,  and  foot  support.  The  test  apparatus  is  shown  in  Figure  1  (adopted  from 
Hammond  and  Roe  1972).  Drivers  were  asked  to  grasp  the  ends  of  the  measurement  rods 
and  to  push  forward  as  far  as  possible.  The  rod  rack  was  moved  laterally  to  span  the 
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space  in  front  of  the  driver.  Testing  was  conducted  with  lap  belt  only  and  with  lap  and 
fixed-length  torso  belts.  The  data  were  analyzed  to  produce  reach  surfaces  and  Figure  2 
depicts  one  such  surface.  Reach  surfaces  are  interpreted  with  respect  to  population 
capability  instead  of  the  capability  of  individuals  with  particular  body  dimensions.  95 
percent  of  drivers  for  a  50/50  male/female  American  driving  population  are  expected  to 
be  able  to  reach  to  push-button  targets  that  are  located  aft  of  the  surface  in  Figure  2 
instead  of  the  reach  capability  of  a  male  driver  who  is  95th  percentile  by  stature. 


Figure  1  SAE  control  reach  measuring  fixture  (Hammond  and  Roe  1972) 


Figure  2  Control  reach  envelop  in  the  vehicle  workspace  (Hammond  and  Roe  1972) 
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Digital  human  figure  models  have  been  used  for  the  driver  workstation  assessments 
(Porter  et  al.  1993,  1995,  and  1998;  Chaffin  2001).  Commercial  digital  human  modeling 
and  simulation  programs  are  commonly  used  to  simulate  reach  to  controls  with  the 
outcomes  of  the  digital  human  based  studies  used  to  assess  control  locations. 

An  experimental  method  was  developed  by  Reed  et  al.  (2003),  which  is  based  on  a 
unified  model  of  reach  difficulty  and  capability.  While  experimented  methods  provide 
insight  into  the  problem,  extrapolating  results  to  populations  that  have  been  used  as 
subjects  remains  problematic.  Furthermore,  extracting  mathematical  models  based  on 
statistical  data  (typically  non-linear  regression)  is  particularly  difficult.  Parkinson  et  al. 
(2003)  compared  the  reach  envelope  determined  entirely  by  the  segment  lengths,  joint 
degrees  of  freedom,  and  joint  ranges  of  motion  to  those  obtained  in  a  laboratory  study  of 
men  and  women.  Accurate  prediction  of  maximum  seated  reach  requires  consideration  of 
balance  and  pelvis  mobility,  neither  of  which  is  closely  linked  to  joint  range  of  motion. 
Sufficient  ranges  of  motion  in  the  shoulder  and  torso  are  also  needed  to  represent 
postures  near  maximum  reach. 

In  our  previous  work  (Abdel-Malek  et  al.  2001a)  we  developed  a  method  for 
understanding  the  workspace  of  upper  extremities  and  Abdel-Malek  et  al.  (2001b) 
developed  a  formulation  for  realistic  posture  prediction.  For  computation  vehicle  design 
systems,  it  is  necessary  to  have  analytical  models  to  verify  driver  reach.  Our  ultimate  goal 
is  to  enable  the  evaluation  of  driver  discomfort  for  different  positions  of  the  controls 
within  our  analytical  description  of  the  reach  envelope. 
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The  understanding  of  trajectory  fonnations  inside  the  driver  reach  is,  to  a  great  extent, 
dependent  upon  the  identification  of  control  barriers  that  exist  as  impediments  to  motion 
and  that  may  hinder  the  execution  of  a  task.  A  rigorous  mathematical  formulation  based 
on  a  kinematic  model  of  the  upper  extremity  will  first  be  introduced.  Because  of  this 
formulation,  we  will  show  that  barriers  inside  the  workspace  are  identified.  More 
importantly,  closed  form  equations  of  the  workspace  will  be  established.  Furthermore,  it 
will  be  shown  that  visualization  of  the  internal  of  the  workspace  provides  a  powerful  tool 
for  barrier  analysis.  The  analysis  includes  restrained  reach  of  drivers,  that  is,  drivers  who 
are  restrained  by  a  non-extending  shoulder  belt  and  the  unrestrained  reach  of  drivers  or 
drivers  restrained  by  a  lap  belt  only. 

Modeling  and  Formulation 

Whereas  the  anatomy  of  limbs  and  their  joints  are  indeed  very  complex  (as  evidenced  by 
the  debate  in  the  literature  on  the  correct  method  for  modeling  joint  motion),  we  will 
employ  a  kinematic  pair  (or  combination  thereof)  as  used  in  the  field  of  robotics  (which 
indicates  a  constrained  kinematics  joint).  For  example,  if  the  resultant  motion  is 
rotational,  the  joint  will  be  modeled  as  a  revolute  joint.  The  effect  of  a  spherical  joint  is 
modeled  as  three  revolute  joints  whose  axes  intersect  at  the  center  of  the  sphere.  Indeed, 
all  anatomical  joints  can  be  modeled  using  basic  kinematic  pairs. 

Using  only  four  parameters  to  describe  one  coordinate  system  with  respect  to  another,  the 
position  and  orientation  of  each  axis  detennine  the  four  parameters  9j,di,ai,aj ,  and 
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hence,  determine  the  resulting  (4  x  4)  transfonnation  matrix.  To  establish  this  matrix,  it 


is  possible  to  observe  that  a  vector  '  v  resolved  in  the  zth  coordinate  system  may  be 


expressed  in  the  (i-l)th  coordinate  system  ('  1  v)  by  performing  four  successive 


transformations  as  follows. 

(a)  A  rotation  about  the  z;1  axis  by  an  angle  of  0l  to  align  the  x;_,  axis  with  the  x;  axis 
(as  shown  in  Figure  3,  x(_j  /  /x;  and  pointing  in  the  same  direction). 

(b)  A  translation  along  the  z;1  by  a  distance  of  di  units  to  make  x;_,  and  x;  aligned. 

(c)  A  translation  along  the  x,  axis  by  a  distance  of  a,  units  to  make  the  two  origins  of  the 
i  and  (/-  1)  systems  coincide  (the  x;  and  the  x;_j  will  also  be  aligned). 

(d)  A  rotation  about  the  x;  axis  by  an  angle  ai  to  coincide  the  two  coordinate  systems. 


Joint i  +  1 


Figure  3  Establishing  coordinate  systems  and  the  four  D-H  parameters 
In  order  to  obtain  a  systematic  representation  of  any  serial  kinematic  chain,  we  define 
q  =  lAi  -  as  the  vector  of  ^-generalized  coordinates  defining  the  motion  of 

a  limb  with  respect  to  another,  where  qj  is  the  individual  DOF  variable.  The  position 
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vector  function  (shown  in  Figure  4)  generated  by  a  point  of  interest  (typically  on  one  of 
the  fingers)  is  written  as  a  multiplication  of  rotation  matrices  and  position  vectors  as 


=Z[fF,R4V  =®(<o  (i) 

«'= i  i= i 

where  both  p  and  R  are  defined  using  the  Denavit-FIartenberg  (D-H)  representation 
method  (Denavit  and  Hartenberg  1955,  Paul  1981,  and  Fu  et  al.  1987)  such  that 


x(q)  = 


x(q) 

p(q) 

z(q) 


cos  qt 

-COS «,  Sing, 

sin  a,,  sing,. 

sin  qj 

COS  a,  COSg,. 

-sin a,  cosqj 

(2a) 

0 

sin  a,. 

cos  a,. 

OS®) 

r  i 

M 

jet 

o 

o 

JB 

sing,.  d,]T 

(2b) 

where  qi  is  the  joint  angle  from  x;_j  axis  to  the  x;  axis,  ci  is  the  shortest  distance 
between  x;._,  and  x;.  axes,  at  is  the  offset  distance  between  z;  and  z(._j  axes,  and  ai  is 
the  offset  angle  from  zi_]  and  z;  axes. 


Figure  4  Definition  of  the  position  vector  function  x(q) 
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The  vector  function  x(q)  characterizes  the  set  of  all  points  touched  by  the  point  of 
interest.  The  aim  is  to  determine  the  envelope  of  this  set.  At  a  specified  position  in  space 
given  by  (xp,yp,zp) ,  Eq.  1  can  be  written  as  a  constraint  function  as 


x{(\)-xp 

y(q)-yP 
*(q  )~zP 


=  o 


(3) 


In  mathematical  terms,  the  expression  defined  by  Eq.  (3)  is  indeed  a  Manifold  with 
boundary  and  cannot  readily  be  visualized,  because  of  the  relatively  large  number  of 
DOFs  used  to  model  the  upper  extremity. 


Joint  limits  (ranges  of  motion)  are  imposed  in  terms  of  inequality  constraints  in  the  form 
of 

q-^q^q-  (4) 

where  q'  and  q[  are  the  lower  and  upper  limits,  respectively,  and  where  i  =  1,  ...  n , 
where  n  is  the  number  of  DOFs.  Note  that  one  joint  could  have  more  than  one  degrees  of 
freedom  (e.g.,  shoulder  joint  is  modeled  as  three  DOF).  In  order  to  include  the  range  of 
motion  in  the  formulation,  we  transfonn  the  inequalities  above  into  equalities  by 
introducing  a  new  set  of  generalized  coordinates  X  =  ...  Ait]T  such  a  joint  variable 

is  parameterized  as 

qt  =<^,i+^f)/2ll(^,C/-^)/2l^nX;.  i  =  l,...,n  (5) 

where  if  sinT,  01,  then  qt  0gj  and  when  sin  0(5U,  then  q;  0gf.  In  order  to  include 
the  effect  of  joint  limits,  it  is  proposed  to  augment  the  constraint  equation  with  the 
parameterized  inequality  constraints  of  Eq.  (5)  such  that 
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G(q)  = 


x(q)~xp 

y(q)-yP 
z(q  )-zp 

q,  -  W;  +  qf  )/2  -  {qf  -  qf  )/2  sin  \ 


=  0,i  =  l,  ...  n 


(6) 


where  q  =[qr  A/]  is  the  vector  of  all  generalized  coordinates.  Note  that  although 

n  -  new  variables  ( )  have  been  added,  n  -  equations  have  also  been  added  to  the 
constraint  vector  function  without  loosing  the  dimensionality  of  the  problem. 


The  Jacobian  (named  after  the  German  Mathematician  Carl  G.  Jacobi)  of  the  constraint 
function  G(q*)  at  a  specific  point  q()  is  the  (3  +  /? )  x  2n  matrix 

G^G/^q"  (7) 

where  the  subscript  denotes  a  derivative.  Note  that  the  Jacobian  is  defined  in 
mathematical  terms  as  the  derivative  of  the  transformation  (Taylor  and  Mann  1972) 
between  x  and  q.  With  the  modified  formulation  including  joint  limits,  the  Jacobian  is 
expanded  as 


G  .  = 

q 


i  jq. 


(8) 


where  qt  =  dq/db  ,  xq  =  dx/dq ,  0  is  a  (3  x  n)  zero  matrix,  I  is  the  identity  matrix,  and 


Xq  = 


-V</,  T?2 


z  z 

<h  <h 


ya 


(9) 


<h  = 


-<(?r-9f)/2)cos4  0 

0  ~{{qu2  -?2)/2)cos^2 

0  0 

0  0 


0 

0 

0 


<{qUn-qLn)l  2)cos2„ 


(10) 
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Because  the  Jacobian  is  not  square  (more  than  three  DOFs),  rank  deficiency  criteria  were 
developed  for  surfaces  that  are  swept  in  space  (Abdel-Malek  and  Yeh  1997).  These 
surfaces  are  called  “barriers”  provide  significant  insight  towards  better  understanding 
reach  and  comfort  within  the  workspace.  Before  addressing  these  criteria,  however,  it  is 
important  to  show  why  the  singularity  of  the  Jacobian  has  a  direct  effect  on  identifying 
barriers  and  human  reach.  A  singularity  (in  the  pure  mathematical  sense)  is  when  the 
Jacobian  has  no  inverse,  i.e.,  a  solution  cannot  be  found.  To  further  explain,  consider  the 
differentiation  of  Eq.  (1)  with  respect  to  time  as 

x  =  xqq  (11) 

where  q  is  the  vector  of  joint  velocities.  Given  the  hand  velocity  (i.e.,  given  x),  the 
calculation  of  q  requires  computing  an  inverse  of  the  Jacobian  x(|.  For  a  singular 

Jacobian,  it  is  not  possible  to  compute  the  required  velocities  for  such  a  path.  It  will  be 
observed  that  such  behavior  is  associated  with  barriers  within  the  reach  envelope  (e.g., 
when  the  arm  is  fully  extended  and  cannot  extend  any  further,  or  when  some  joints  in  the 
ann  have  reached  their  limits). 

We  will  use  the  idea  of  a  singular  Jacobian  to  identify  all  barriers  inside  and  on  the 
boundary  of  the  workspace.  Because  the  Jacobian  is  nonsquare,  we  define  such  barriers 
as  a  subset  of  the  workspace  at  which  the  Jacobian  of  the  constraint  function  of  Eq.  (7)  is 
row  rank  deficient;  i.e.,  banders  defined  by  a  subset  of  the  reach  envelope  iW  and 
characterized  by 

dW  c|RankG(|,(q*)<(«-2),  for  some  q  with  G(q  )  =  oj  (12) 
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where  k  is  at  least  (3  +  ri) .  Because  of  the  form  of  the  Jacobian  characterized  by  Eq.  (12), 


three  distinct  conditions  arise: 

(1)  Type  I  singularity  sets:  If  no  joints  have  reached  their  limits,  the  diagonal  sub-matrix 
qt  is  full  row  rank.  Therefore,  the  only  possibility  for  G^,  to  be  row-rank  deficient  is 

when  the  block  matrix  xq  is  row  rank  deficient.  Type  I  singularity  set  is  defined  as 

S{1)  =  jp  e  q  :  Rank[xq]  <  3,  for  some  constant  subset  of  qj  (13) 

where  p  is  within  the  specified  joint  limit  constraints  and  may  contain  joints  that  are 
functions  of  others  or  constant  values. 

(2)  Type  II  singularity  sets:  When  certain  joints  reach  their  limits,  e.g., 

Sq1""  =  [^ihmit,^-mit,^I"llt]r ,  the  corresponding  diagonal  elements  in  the  matrix  qA  will  be 

equal  to  zero.  Therefore,  the  corresponding  matrix  is  subjected  to  the  rank-deficiency 
criterion,  where  Gq,  will  take  on  the  following  form 

\  -  x*  \  xft  -  x? 

0  ...  1  0  0  ...  0 

q*  0  ...  0  1  0  ...  0 
0  ...  0  0  1  ...  0 

and  where  the  three  columns  pertaining  to  \q ,  \q  ,  and  x?  are  removed  such  that  the 
rank  deficiency  criteria  are  applied  again.  From  the  foregoing  observation,  the  second 
type  of  singular  sets  are  formulated.  Define  a  new  vector  ^q 1111111  =  ^/'"'"‘.c/'"1"1,^)"1111  j  , 

which  is  a  sub-vector  of  q  where  1  <  dirnGl"mitll(n-3). 

The  type  II  singularity  set  is  defined  as 

Si2)  =  jp  =  [p  u^qlimit] :  Rank[xq(w,<3qlim“)]  <  3,  for  some  peq,  dim(<3qhmit)  <  (n  -3)| 
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(15) 

where  p  is  the  singular  set  as  a  result  of  applying  the  rank  deficiency  criteria  to  Eq.  (14). 
(3)  Type  III  singularity  sets:  are  all  sets  that  are  composed  of  the  combination  of  joints 
at  their  limits  and  is  defined  by: 

S(3)  p  ©<T|limit  H[g,limit,^Mt,...]S;  where  i  © j  (16) 

Barriers  are  identified  by  substituting  the  sets  p;  characterized  by  Eqs.  (13,  15,  and  16) 
into  the  accessible  set  x(q) ,  which  yields  the  equation  of  a  surface  that  can  be  readily 
shown.  This  surface  is  indeed  a  barrier  associated  with  a  generalized  variable  that  has 
reached  its  limit.  Determining  joint  angles  of  the  upper  extremity  given  a  specific 
position  and  orientation  is  usually  defined  as  the  inverse  kinematics  problem  in  the 
robotics  literature  (Fu,  et  al.  1987).  Motion  from  one  configuration  to  another  along  a 
trajectory  sometimes  requires  halting  the  motion  and  changing  the  inverse  kinematics  in 
order  to  proceed  with  the  motion.  An  example  of  this  occurs  when  attempting  to  reach  a 
point  located  behind  one’s  shoulder.  Starting  with  one  trajectory  may  become  very 
uncomfortable  because  of  joint  limits,  while  trying  another  trajectory  becomes  simpler 
and  uninterrupted.  Similarly,  reaching  a  doorknob  and  turning  sometimes  is  difficult  to 
complete  and  requires  orienting  the  initial  hand  configuration  in  a  different  posture. 

These  barriers  due  to  singular  sets  identified  by  Eqs.  (13,  15,  16)  may  admit  motion  only 
in  one  nonnal  direction,  and  hence  are  called  impediments  to  motion. 

Motion  on  a  Barrier 

To  better  understand  when  the  hand  may  or  may  not  cross  barriers  under  given 
conditions,  we  explore  the  barrier’s  kinematic  properties.  We  propose  a  criterion  that  is 
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based  on  nonnal  acceleration  at  a  point  on  a  barrier,  such  that  crossability  is  achieved  if 
the  barrier  admits  a  nonnal  acceleration  in  one  direction  or  another.  A  point  on  a  barrier 
admits  motion  normal  to  the  surface  in  either  direction  depending  on  the  difference  in 
acceleration  components  (defined  by  the  indicatory),  such  that 

2 


V  =  an 


(17) 


where  vt  is  the  tangential  velocity,  an  is  normal  acceleration,  and  1 / p0  is  the  normal 
curvature  of  the  barrier  with  respect  to  the  tangent  direction  of  vt  ( po  is  the  radius  of 
curvature).  The  need  for  formulating  the  problem  in  terms  of  velocities  and  accelerations 
will  become  apparent,  as  the  resulting  expression  for  the  indicator  //  will  be  independent 
of  acceleration  values,  but  will  be  a  quadratic  fonn  that  has  definiteness  properties.  A 
point  on  a  singular  surface  will  have  no  acceleration  if  the  quantity  rj  computes  to  null. 


For  a  singular  parametric  entity  f(,)(u<!))  e  R3  (where  u  is  a  vector  representing  the 
remaining  joint  variables — those  not  constant),  and  in  the  field  of  differential  geometry, 
the  First  Fundamental  Form  (Farin  1993)  is  denoted  by  I  ,  where 

u<0  =  [u  vf  =  [y(.  J  ,  and  is  defined  as 

Ip^rf„rf>  (18) 

where  fu  =  ff/fu  .  The  Second  Fundamental  Form  is  defined  as 

IT  =  (^ir[Nrfl  fiu  (19) 

P  L  -luu 

or  expanded  to 

II/;  =  Nrf Ulldu2  +  2Nrf Uvdu  dv  +  N%,dv2  (20) 
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where  N  is  the  vector  normal  to  the  singular  surface  and  fMV  =  rff /chick .  The  Nonnal 
Curvature  Ko  of  a  parametric  surface  at  a  configuration  qo ,  in  the  direction  of  du/dv , 
can  then  be  defined  as  the  ratio  (Farin  1993) 


O  T 

Po  1P 


(21) 


In  order  to  detennine  the  kinematics  quantities,  we  define  the  Time-Modified  First  and 
Second  Fundamental  Forms  as 


i;  =  urf„rfnu  (22) 

(23) 


such  that  the  normal  curvature  can  still  be  defined  as 


Po  IP  K 


(24) 


For  a  singular  surface  f  *' 1  (u(,) ) ,  the  derivative  using  the  chain  rule  is  fuii .  Similarly,  for 
the  general  singularity  <t>(q) ,  the  derivative  is  Oqq  .  Therefore,  at  an  instant  of  time,  the 

tangential  velocity  in  terms  of  f  or  <1>  at  any  point  on  the  barrier  is 

v,=fuu  =  Oqq  (25) 


If  joint  limits  are  considered,  then  the  derivatives  can  be  written  as  u  =  and  q  =  q;i . 
The  squared  norm  of  the  velocity  is 


V  =v>(  =urf«rfuu 


(26) 


which  is  equal  to  the  Time-Modified  First  Fundamental  Form  1^  of  Eq.  (22).  Therefore, 
I'  can  be  written  as 


Substituting  l/po  into  r/  yields 


(27) 


d  =  «„- v, 


K 

i' 


=  a- 


II' 


(28) 
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When  we  consider  Jacobian  singular  surfaces  an  =  Nrx  =  qr  j^NrO  J  q  .  Since  an  is  in 

tenns  of  q  and  to  express  II in  terms  of  q ,  it  was  shown  that  the  velocity  vector  on 
singular  surface  can  be  written  as 


where  E  = 


1 

0 


*i  =  [Ef„]  1  EOqq  (29) 

0  0 
1  0 


if  the  first  and  second  rows  of  fu  are  independent;  E  = 


1  0  0 
0  0  1 


if  the  first  and  third  rows  of  fu  are  independent;  E  = 


0  1 
0  0 


0 

1 


if  the  second  and  third 


rows  of  fu  are  independent.  Therefore  the  crossability  criteria  can  be  expanded  into  a 
quadratic  form  by 

V  =  an-ll'p=  qQq  (3°) 

where 

Q  =  rNrol  -O^bT  Nrf  1  BO  (31) 

L  Jqq  4  L  Juu  1 

and  B  is  the  generalized  inverse  of  fu  defined  by 

B  =  [Ef„]  1 E  (32) 

The  matrix  Q  is  calculated  by  substituting  q0  into  Eq.(31).  The  criteria  is  that  if  Q  is 

indefinite,  the  interest  point  of  the  driver  can  admit  normal  movements  along  either 
direction  of  N  .  The  surface  is  crossable.  If  Q  is  semi-definite  the  surface  is  non- 
crossable. 


When  we  consider  the  joint  limits  the  component  of  the  normal  acceleration  is  then 
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an  =  Nrx  = 


(33) 


where 


H'=tf[NHA+£ 

i= 1 


^/(NrO) 

dq, 


kk 


(34) 


and  the  quadratic  form  is  written  in  tenns  of  the  Q  matrix  as 

Q*  =  H'-q>[Br[Nrf]jiB®1ql 


(35) 


If  Q  is  indefinite  the  surface  will  definitely  be  crossable.  However  when  Q  is  either 
positive  semi-definite  or  negative  semi-definite,  the  singular  surface/curve  may  still  be 
crossable.  To  address  this  case,  we  propose  the  projection  of  a  variational  movement 
Ox  =  x  Sqf  due  to  5qj  onto  the  normal  direction  N  such  that  the  normal  component 

o"  =  N  T\5q,  (36) 

determines  admissible  normal  movement,  where 

_  +1  if  q,  is  at  lower  bound 

oq,  =  \  (37) 

[-1  /  qi  is  at  upper  bound 

Positive  value  of  a  in  Eq.  (37)  indicate  the  tip  of  the  middle  finger  can  admit  movement 
in  the  positive  direction  of  N  . 


Another  situation  arises  when  the  nonnal  vector  N  is  perpendicular  to  x?  .  In  this  case, 
a  evaluates  to  zero.  Therefore  we  define  the  curvature  difference 


fi<7,  [ Nrx~|  Sqt  fiurrNrfl  fiu 

L  Mat _ L  Juu 

Sq\lqxqSqi  ^ufurfufiu 


(38) 


We  have 
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Ju  =  BxqSqi 


(39) 


Plugging  Eq.(39)  into  Eq.(38)  yields 

<J?)[Nrx]  Sq,-Sq,xTt Br  [Nrf  ]  B xjq, 

K.-K,= - 1 - iMi - - — L — ii! — : -  (40) 

dq^qXqSq, 

The  denominator  of  the  right-hand  side  of  Eq.(40)  is  always  positive,  the  sign  is 
dependent  on  the  numerator.  The  numerator  can  be  reduced  to 

<u  =  KSqf  (41) 


where 

^  =  x^BrrNrxl  Bx  -Tr^xl  (42) 

"  L  Auu  "i  L  J 

If  K  >  0 ,  the  tip  of  the  middle  finger  can  admit  movement  into  the  positive  direction  of 
N  .  If  K  <  0 ,  the  tip  of  the  middle  finger  can  admit  movement  into  the  negative  direction 
of  N  .  When  Q  is  either  positive  semi-definite  or  negative  semi-definite.  The  final 
criteria  will  be: 

(1) .  If  N  is  not  perpendicular  to  xq  ,  then  a  in  Eq.(36)  must  be  evaluated. 

(2) .  If  N  is  perpendicular  to  xq  ,  then  K  in  Eq.(42)  must  be  evaluated. 

If  any  of  a  or  K  has  a  different  sign  than  the  nonzero  eigenvalues  of  Q  ,  the  singular 
surface/curve  is  crossable;  If  a  and  K  have  the  same  sign  as  nonzero  eigenvalues  of 
Q  ,  the  singular  surface/curve  is  non-crossable. 
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The  Analysis  of  Restrained  Drivers 

Finding  the  Reach  of  Drivers 


Consider  the  shoulder  and  forearm  modeled  as  a  4-DOF  system,  where  the  spherical  joint 
at  the  shoulder  is  modeled  as  three  intersecting  revolute  joints  and  the  elbow  as  a  revolute 
joint.  This  is  consistent  with  published  results  except  that  we  have  not  considered  the 
wrist  joint  (which  is  an  additional  3  revolute  joints)  and  we  have  limited  the  motion  of  the 
glenohumeral  joint  to  spherical.  It  should  be  noted  that  this  model  does  not  include  the 
scapulothoraic  motion  of  the  shoulder  joint  (i.e.,  the  additional  two  translational  DOFs  of 
the  scapulothoraic  are  not  taken  into  consideration).  Furthermore,  we  have  modeled  the 
spherical  joint  as  three  revolute  joints  intersecting  at  one  point,  a  practice  commonly 
made  in  modeling  to  represent  spherical  joints.  It  should  be  noted  that  the  most  difficult 
and  the  least  successful  modeling  of  a  major  articulating  human  joint  has  been  the 
shoulder  because  of  the  difficulty  in  modeling  the  complex  kinematic  connectivities  as 
well  as  the  anatomical  complexity  of  the  region.  Figure  5  depicts  the  joint  motions  to  be 
modeled  where  each  joint  is  given  an  independent  coordinate  qt  where  the  equivalent 
kinematic  skeleton  of  the  system  is  depicted  with  the  z-axis  located  per  the  D-H 
representation  method. 
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Figure  5  Shoulder  and  arm  and  the  corresponding  DOF 


In  the  following  analysis,  a  point  on  the  tip  of  the  middle  finger  as  shown  in  Figure  6  will 
be  tracked.  The  dimensions  of  the  arm  are  different  for  different  percentile  humans.  For 
95th  percentile  men  the  dimensions  areZl  =  40.64cw  ,  L2  =  30.48 cm  ,  L3  =  20.32 cm  , 
nevertheless,  our  method  allows  for  any  anthropometric  data  to  be  used  in  the 
formulation  and  is  not  percentile  dependent.  In  the  field  of  kinematics,  the  motion  of  a 
spherical  joint  with  three  DOFs  can  be  modeled  as  three  independent  revolute  joints 
having  their  axes  intersecting  at  a  single  point  as  shown  in  Figure  6.  Note  that  the  point 
on  the  middle  finger  is  shown  located  at  the  position  4v  =  [0  0  50.8]r  as  resolved  in 
the  fourth  coordinate  frame. 
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Figure  6  Kinematic  modeling  of  the  restrained  driver 


Consider  the  following  joint  limits  imposed  on  the  model  of  Fig.  6  for  a  95th  percentile 
male:  050°  Oq\ 090,  0510°  Oq2  0120°,  05O°  0^3  090°,  and  -150°  <  q4  <  0°.  Using 
the  Denavit-Hartenberg  representation  method,  the  position  of  the  point  of  interest  is 
given  by  Eq.  (1)  as 


X(q) 


(40.64  + 50. 8cos <y4 ) sin  ql  sin qi  +cos qx(cos ^(-40.64- 50. 8cos gjsin q2  -50.8cos q2  sin t/4) 
cosg,3(-40.64-50.8cos ^4) sin q1  sing,  +cos g[(-40. 64-50. 8cos g4)sin q2  -50.8cos q2  sin g,  sin g4 
cos q2  cosg3(40. 64  +  50. 8cosg4)-50.8sin q2  sing4 


(29) 
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Inequality  constraints  are  parameterized  as  ql  =  n  /  2  sin \ ,  q2=  ft  /  36  +  Tin  /  36  sin  /l, , 


<73  =  ;r  /  2 sin  /l3 ,  q4  =  -5;r / 1 2  +  5zr  / 1 2  sin  /.4 .  The  rank  deficiency  criteria  applied  to  the 
resulting  (3  <M)  Jacobian  matrix  yields  44  singular  sets  that  are  listed  in  Appendix  A. 
Those  singular  sets  are  now  substituted  into  Eq.  (29),  whereby  parametric  surface 
equations  are  readily  plotted  and  visualized.  For  example  several  cross  sections  of  the 
final  driver  reach  envelope  are  shown  in  Figure  7. 


Figure  7.  The  cross  sections  of  the  restrained  driver  reach 
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Crossabilitv  Analysis  of  Singular  Surfaces 


A  planer  y0oz0  cross  section  view  of  the  reach  envelope  at  ( x0  =  0  )  is  shown  in  Figure  8, 
where  the  intersection  of  the  singular  surfaces  with  plane  traces  curves  that  are  called  the 
singular  curves  and  are  denoted  by  VFJ ,  i  - 1,...,44 . 


Figure  8  Cross  Section  of  Singular  surfaces 

To  demonstrate  the  crossability  analysis,  consider  a  point  p,  on  surface  y/n  ,  which  has 
q  =  [0  Inti  nil  -7vllf  and  k  =  [0  n/2  nil  -arcsin(l/5)]r .  Evaluating  the 

nonnal  N  from  the  basis  of  the  null  space  of 

63.8372  -39.8982  0" 

0  0  0 

0  0  0 

0  -65.1535  0 

at  Pj  yields  the  normal  vector  N  =  [0  0  I  ]7.  Evaluating  the  matrix  H  of  Eq.  (34) 
yields 


r  n7 

_®qq>_  = 
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The  singular  surface  is 


0 

0 

0 

0 


0 

50.9811 

0 

0 


0  0 

0  0 

-31.9186  0 

0  -72.3674 


f(u)  =  [(40.64  +  50.8 cos qA) sing,  +25.4cos(/,  sing4  cosg,(-40. 64-50. 8 cos (/4)  + 25 .4 sing,  siru/4  -43 .9941  sing,,] 7 


and  its  Jacobian  matrix  can  be  defined  by 


fu(0,-7T/2)  = 


40.64  0 

-25.4  -50.8 
0  0 


Therefore  the  matrix  E  is 


E  = 


1  0  0 
0  1  0 


and 


B  =  [Ef„]  1 E  = 


0.0246063 


0 


0 


-0.0123031  -0.019685  0 


The  matrix  Q  of  the  quadratic  form  is  evaluated  as 

"0  0  0  0 


Q  = 


0  50.9811  0  0 

0  0  -31.9186  0 

0  0  0  0 


The  eigenvalues  of  Q*  are  evaluated  as  {50.981 1  -31.9186  0  0}  ,  which  indicates  an 


indefinite  quadratic  fonn.  Therefore  it  is  a  crossable  surface  at  point  Pj .  Crossability  at 
this  point  means  that  a  driver  with  a  posture  characterized  by  the  singular  set  su ,  would 

allow  the  hand  to  cross  that  surface  (represented  by  a  singular  curve).  This  is  the  first  step 
towards  better  understanding  comfort,  dexterity,  effort,  and  energy  zones  associated  with 
driver  reach. 
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To  further  demonstrate  this,  consider  a  point  p2  on  surface  i//,  5 ,  which  has  a 
configuration  of  q  =  [-;r/  2  0  0  0]r  and  k  =  \-n!2  arcsin(-l/23)  0  ;z72]r.The 

nonnal  vector  at  this  point  is  calculated  as  N  =  [0  0  l]r  .  The  matrix  Q  is  computed 
and  its  nonzero  eigenvalue  is  -22.5778  .  Since  Q  is  negative  semi-definite,  the  singular 
surface  y/X5  is  non-crossable. 

Consider  a  point  p3  on  surface  i//4 ,  which  is  a  Jacobian  singular  surface  and  has  a 
configuration  q  =  [^ / 2  2n!2>  0  -;z7  2]r  and  k  =  \n!2  nil  0  -arcsin(l/5)]r . 

The  normal  vector  is  N  =  [0  -0.5  0.866025]7  .  The  nonzero  eigenvalue  of  Q  is 

obtained  as  81.5697.  Since  the  matrix  is  positive  semi-definite,  the  supplementary  criteria 
must  be  computed  to  ascertain  crossability.  For  this  surface,  the  joint  variables  q]  and  q2 

are  at  their  limits.  The  variational  nonnal  movement  of  <7n  for  Sq]  and  Sq2  are  also 

calculated  as  0  and  -40.64  respectively.  For  5qx  the  value  of  an  is  0  and  we  can 

calculate  the  Eq.(42)  and  K  =  -30.2976 .  The  signs  of  on  for  Sq2  and  K  for  8qx  are 

different  from  the  sign  of  nonzero  eigenvalue.  Therefore,  the  surface  i//4  at  this  point  is 

crossable.  Use  the  same  procedure  to  check  the  points  on  other  surfaces  and  the  final 
admissible  normal  movement  directions  are  shown  in  Figure  9. 
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:  Crossable  Surfaces 


:  Movable  Directions 


Figure  9  Crossable  and  non-crossable  singular  curves 

The  Unrestrained  Reach  of  Drivers 

Consider  the  upper  body  modeled  as  a  7-DOF  system,  where  the  upper  extremity  is  the 
same  as  the  restrained  model  in  Figure  6.  In  order  to  account  for  the  upper  torso  motion 
for  an  unrestrained  driver,  we  propose  the  addition  of  3  degrees  of  freedom  representing 
motion  at  the  waist.  The  torso  part  has  three  DOFs  (spherical  joints)  and  they  are 
intersected  into  one  point  characterizing  a  seated  driver  with  spherical  articulation  at  the 
waist.  Consider  the  following  joint  limits  imposed  on  the  model  of  Figure  10  for  a  95th 
percentile  male  model:  -45°  <  ql  <  45 0 ,  0°  <q2<  30° ,  -30°  <  q3  <  30° ,  -90°  <  q4  <  90° , 

-110°  <  q5  <  120° ,  -90°  <q6<  90 0 ,  and  -150°  <q7<  0° . 
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Figure  10.  Kinematic  modeling  of  the  unrestrained  driver 


According  to  the  D-H  method,  the  tracked  point  of  interest  is 

x(q)[l]  =  25 c,c3s2  +  25s,s3  +  10c3s,  -lOq^s,  +  \6c6(c5(cisx  -cls2si)-(c4(clc2s2  +  + 

+  16(CjC,C4  +  (c,C35,  +5153)s4S6  +  20(c7(c6(c,(c35j  -C1S2S3)-(c4(c1C3,5'2  +51S3)-C,C254)s5)  + 
(c,C2C4  +  [CyC2S2  +  X^5'3)54)Xg)  +  (  C3 (c4 (C,C352  +  SyS2 )  C^C2S4 )  (c3^  CyS2S2 )^7 ) 


x(q)[2]  =  25c3s1.s'2  +25xjX3  -q,s3  +10(-c,c3  -x,x2x3)  +  16c6(c5(-c1c3  -x1x2x3)-(c4(c3x,x2  -c,x3) 


-^54)55) +  16(c2c4x,  +(c3XjX2  -CjX3x4)x6  +  20(c7(c6(c5(-c,c3  -x1x2x3)-(c4(c3x,x2 

-c1x3)-c2x,x4)x5)  +  (c2c4x1  +(c3x,x2  -c1x3s4)x6)  +  (-c5(c4(c3x1x2  -c1x3)-c2x1x4)- 

(-CiC3  -.SVS'2.S'3),S’5).S'7) 
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x(q)[3]  =  25 c2c3  -10c2s3  +16 c6(-c2c5s3  -(c2c3c4  +  s2s4s5 )  + 1 6(-c4s2  +c2c3s4)s6  +20(c7(c6(-c2c553 


-  (c2C3C4  +  J2,VS  )  +  (— C4s2  +  C2C354  )s6  )  +  (— Cj  (C2C3C4  +  52^4  )  +  C2S3S5  )s7  ) 


Figure  1 1  shows  the  cross  section  of  the  unrestrained  driver  reach  in  the  car.  A  x0z0  cross 

sectional  plane  showing  the  intersections  of  the  singular  surfaces  with  the  plane  as 
singular  curves  is  shown  in  Figure  12. 


ACCELERATOR 

heel  point 


Figure  11.  The  cross  section  of  the  unrestrained  driver  reach  ( y0  =  -10  ) 
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Figure  12.  The  cross  section  at  x0  =  0 

We  follow  the  procedure  for  determining  the  crossability  of  singular  surfaces  and  obtain 
the  crossabilities  of  the  boundary  surfaces  of  the  unrestrained  driver  reach  in  Figure  13. 


Figure  13.  Boundary  surfaces  and  their  non-crossabilities  (x0  =  0) 
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Obviously,  unrestrained  driver  reach  gives  rise  to  significantly  more  redundancy  as  there 
are  significantly  more  singular  curves.  It  is  also  observed  that  far  more  crossable  surfaces 
in  both  directions  are  observed  in  the  unrestrained  case,  as  expected,  because  of  the 
redundant  7  DOF  model. 

Conclusions 

A  formulation  for  visualizing  restrained  and  unrestrained  driver  reach  envelope  has  been 
presented,  particularly  the  identification  of  singular  surfaces  within  upon  which  Jacobian 
singularities  occur.  The  work  is  aimed  better  understanding  human-machine  interaction 
from  a  human  perspective,  where  the  ultimate  goal  is  to  delineate  regions  and  zones  of 
common  and  quantifiable  characteristics  such  as  comfort,  effort,  energy,  etc.  This  initial 
stage  has  focused  on  the  understanding  of  barriers  within  the  envelope  upon  which 
specific  postures  may  or  may  not  allow  for  crossability  from  kinematic  point  of  view.  It 
was  shown  that  restrained  drivers  whose  upper  extremity  is  modeled  as  a  4DOF 
kinematic  linkage  has  far  more  non-crossable  (i.e.  crossable  surfaces  in  one  direction  at  a 
particular  singular  configuration)  than  unrestrained  drivers  modeled  as  a  7DOF  kinematic 
linkage.  It  was  also  shown  that  crossability  at  a  singular  curve  within  the  reach  envelope 
can  now  be  determined  in  closed  form  and  based  as  a  first  but  necessary  step  towards 
delineating  zones  and  regions  within  the  envelope  associated  with  driver  quantifiable 
properties. 
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Appendix  A 

si  =  {qx  =-71/2, q2  =  — 1  lzr / 1 8} ,  s2  =  {qx  =  ;r/2, q2  =  —  1  lzr / 1 8} , 

V  =  {qx  =-K/2,q2  =  2zr / 3} ,  54  ={qx  =njl,q2  =  2zr/3}  , 

V  ={<7i  =-tf/2,03  =-7tI2}  ,s6={qx  =  -n/2,q3  =n!2)  , 
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s7  =  {qi  =  x/2,q3  =  — tt  /  2}  ,jg  ={qx  =  x/2,q3  =  x / 2} , 

sg  ={q2  =-\\n 1 18,^  =  7T / 2}  ,s10  ={q2  =  -lbr/18,<73  =  — tt / 2} , 

Sji  =  {^2  =  2x  13,q2  =  7T / 2}  ,512  =  {g2  =  2tt / 3 ,q3  =  -x  I 2} 

•5,3  ={^i  =  -x/2,q4=-5x/6},  s14  =  {^  =  x/2,q4  = -5x / 6} , 

*i5  =  {<h  =-K/2,q4  =0},  j16  ={^  =  x/2,q4=0), 
sv  ={q2  =  -1  l;r/18,g4  =  -5^/ 6}  ,  slg  =  {g2  =  -1  Itt / 1 8, q4  =  0} , 
s19  =  {<?2  =2^/3 ,g4  =  -5;r / 6}  , s20  =  {<?2  =2x13, q4  =0}, 
s21  ={^3  =~xl2,q4  =-3x16),  s21  =  {q3  =x/2 ,q2  =-5x16}  , 

*23  ={^3  =-x/2,q4  =0},  s24  ={^3  =nl2,q4  =0}, 

525  =  =  -x/2 ,q4  =-^rccos(-^-)|,  s26  =  =  x/2,q4  =  -^rccos(-^)| , 

*27  =  {^2  =  “1  1^/18, =  0}  ,s28  =  { q2  =  2x13, q3  =  0}  , 

f  4  1  f  4  1 

*29  =  \q2  =  2x  13,  q4  =  -Arc  cos(-— )  m30  =  j  q2  =  -1  br/l  8,  q4  =  -Arc  cos(-— )  >, 

*31  =  {^3  =  -x  1 2, q2  =  -x  1 2}  ,s32  =  }q3  =  -x  / 2, q2  =  x  1 2}  , 

*33  =  {<?3  =  x !2,q2  =  —x 1 2} ,  s34  ={q2  =  x !2,q2  =  x  1 2} , 

*35  ={<h  =0,q4=-5x/6},  s36  ={q2  =0 ,q4  =-5x/6), 

*37  =  {<?2  =^/2,q4=-5x  /  6} ,  s3i={q2  =-x/2,q4=-5x  /  6} , 

*39  ~  {q2  =x/2,q4=0}  ,s40  =  {q2  =0,q4=0}, 
s4l={q2  =  x/2 ,q4  =  0} ,  s42  =  [q2  =  -x/2,q4  =  0} , 

*43  ={^3  =x!2,q4  =0},*44  ={q3  =0,q4  =0} 
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